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ABSTRACT 

TypeT X-ray bursts are thermonuclear flashes that take place on the surface of accreting neutron stars. The wait time between 
consecutive bursts is set by the time required to accumulate the fuel needed to trigger a new burst; this is at least one hour. 
Sometimes secondary bursts are observed, approximately 10 min after the main burst. These short wait-time bursts are not yet 
understood. We observed the low-mass X-ray binary and X-ray burster EXO 0748-676 with XMM-Newton for 158 h, during 
7 uninterrupted observations lasting up to 30 h each. We detect 76 X-ray bursts. Most remarkably, 15 of these bursts occur 
in burst triplets, with wait times of 12 min between the three components of the triplet, Tl, T2, and T3. We also detect 14 
doublets with similar wait times between the two components of the doublet, Dl and D2. We characterize this behavior to try 
and obtain a better understanding of bursts with short wait times. We measure the burst peak flux, fiuence, wait time and time 
profile, and study correlations between these parameters and with the persistent flux representing the mass accretion rate, (i) 
For all bursts with a long wait time, the fiuence is tightly correlated with the wait time, whereas burst with short wait times 
generally have higher fluences than expected from this relationship; (ii) Wait times tend to be longer after doublets and triplets; 
(iii) The time profile of single bursts, SI, and of the first burst in a double or triple burst, Dl and Tl, always contains a slow 
component which is generally absent in the D2, T2 and T3 bursts; (iv) The peak flux is highest for SI, Dl and Tl bursts, but 
this is still a factor of 7 lower than the highest peak flux ever seen for a burst in this system; (v) The persistent flux, representing 
the mass accretion rate onto the neutron star, is about 1% of Eddington, which is among the lowest value so far measured 
for this system. The amount of energy per gram of accreted mass liberated during bursts is consistent with a fuel mixture of 
hydrogen-rich material. The characteristics of the bursts indicate that possibly all bursts in this system are hydrogen-ignited, in 
contrast with most other frequent X-ray bursters in which bursts are helium-ignited, but consistent with the low mass accretion 
rate in EXO 0748-676. Possibly the hydrogen ignition is the determining factor for the occurrence of short wait-time bursts. For 
example the 12 min wait time may be associated with a nuclear beta decay timescale. 
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1. Introduction 

Type I X-ray bursts are due to unstable burning of hydro- 
gen and/or helium on t he surface on an ac creting neutron 



star (for a review see iLewin et al. 1993)- Fuel accumu- 



lates for hours to days, and then ignites as soon as the 
pressure and temperature conditions for thermonuclear re- 
actions are reached. Since the nuclear reaction rates de- 
pend strongly on the temperature, the ignition leads to 



a runaway process, and the fuel burns explosively in a 
rapid ~10-100 s burst releasing typically 10 38-39 ergs. 
Two bursts from a given source are usually separated by a 
few hours. However, doublets with burst intervals as short 
10 minutes have been observed from several sources 



a.s 
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including EXO 0748-676 (Gottwald et al. 


19861). GS 0836- 


429 (|Aoki et al. 1992ri. 4U 1608-522 ( 


Murakami et al. 


1980b 4U 1636-536 Johashi et al. 1982 


; IPedersen et al. 


1981. 4U 1705-440 (lLanemeier et al.lll987) and XB 1745- 


248 (Inoue et al. 1984f). A burst interval as short as 50 s 



to: L. Boirin, e-mail: has been seen in MXB 1659-29 (jWijnands et all 120021 ) 
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Table 1. XMM-Newton observations of EXO 0748-676 performed between September 19 and November 12, 2003. We 
indicate the observation identification and revolution numbers, the EPIC PN start and exposure times, the number of 
X-ray bursts and bursting events observed (for singlets, doublets and triplets separately and in total), and the mean 
event wait time. 



Obs-ID 


Rev. 


Start (UT) 
Month Day h:m 


^cxp 


Bursts 


Singlets 


Doublets 


Triplets 


Events 


Mean event wait time 
(h) 


0160760101 


692 


Sept. 


19 


13:37 


24.6 


10 


6 


2 





8 


3.20 


0160760201 


693 


Sept. 


21 


13:38 


25.1 


14 


6 


1 


2 


9 


2.65 


0160760301 


694 


Sept. 


23 


10:42 


30.0 


14 


7 


2 


1 


10 


3.02 


0160760401 


695 


Sept. 


25 


17:29 


20.4 


a g 


3 


a 3 





a 6 


3.49 


0160760601 


708 


Oct. 


21 


10:02 


15.2 


8 


3 


1 


1 


5 


2.99 


0160760801 


710 


Oct. 


25 


19:19 


17.3 


9 


3 


3 





6 


2.71 


0160761301 


719 


Nov. 


12 


08:24 


25.2 


12 


5 


2 


1 


8 


3.19 






Total 




157.9 


76 


33 


14 


5 


52 


Mean 3.02 



a including a burst detected by RGS only, a few seconds before EPIC cameras were turned on. 



Lewin et al. ( 19761 ) report SAS-3 observations of three con- 
secutive bursts from MXB 1743-28 with burst intervals of 
18 and 4 minutes. However, these are observations of a 
field with several known bursters, and therefore there is a 
slight chance that these three bursts may not be attributed 
to one source. 

Bursts separated by short time intervals are not well 
understood in the context of the classical thermonuclear 
flash model. Since the interval is too short to accrete 
sufficient material to fuel a thermonuclear burst, sec- 
ondary bursts are thought to be due to burning of resid- 
ual fuel that did not burn during the primary burst. 
This involves incomplete nuclear burning, fuel storage and 
a mixing mechanism with the freshly accreted material 



( Fuiimoto et~allll987[ ) 



EXO 0748-676 is one of the sources where the dou- 
ble burst phenomenon was best studied thanks to sev- 
eral EXOSAT observations carried out at different per- 
sistent fluxes of the system. In 1985, 26 bursts were de- 
tected, including four doubl ets with a burs t sepa ration 



of the order of 10-20 min (jGottwald et all Il986l ). The 



doublets only occurred when the persistent flux was low. 
As the persistent flux increased, the doublet phenomenon 
stopped, the wait time of the single bursts increased, 
and their shape changed from a "slow" (long tail) to a 
"fast" profile sometim es showing photospheric expansion. 
Gottwald et al.l {l986) hypothesized that these variations 
could be caused by the flashes changing from hydrogen- 
triggered hydrogen-helium flashes at low accretion rates to 
helium-dominated flashes at high accreti on rates, a sce- 
nario previously theoretically studie d by iFuiimoto et al 



( 198lh (see also lAvasli fe Joss 1982 ) who outlined three 
ways to trigger shell flash es depending on the accretion 
rate. lGottwald et alJ (|l986l ) speculated that double bursts 
could be a feature of hydrogen-triggered hydrogen-helium 
flashes. 

In another EXOSAT observation in 1986, 11 bursts 
were observed and showed a regular pattern with a long 
recurrence time always followed by a short one, reminis- 
cent of the double burst phenomenon, but with longer sep- 



arations, in the range 20-70 min (jGottwald et al.l 119871 ). 
The wait time to a burst and the total emitted energy in 
that burst displayed a linear relation but with an offset 
energy at zero burst interval that was interpreted as in- 
complete consumption of fuel in the primary burst and 
its subsequent consumption in the secondary burst. The 
amount of unburned fuel was estimated to be 10-15% of 
the total available nuclear energy. 



Cottam et al.l (|2002[ ) reported redshiftcd spectral lines 



identified with O and Fe transitions during X-ray bursts 
from EXO 0748-676 observed with RGS on XMM-Newton, 
implying a gravitational redshift z = 0.35. Combining this 
redshift measurement wit h the peak flux of t he radius ex- 
pansion burst observed by I Wolff et a l. (2005) and the flux 
and color temperature observed dur ing the coo ling tail 
of bursts from EXOSAT and RXTE, lOzell (|2006l) derived 
limits on EXO 0748-676 compact object mass and radius 
which ruled out several soft equations of state for the neu- 
tro n star interio r. 

o3 (2006) further derived a lower limit to the 
EXO 0748-676 distance of 9.2 ± 1.0 kpc (see note ED in 
Sect. [6]). The source distance was previously estimated 
from bursts showing ph otospheric expansion to ra nge be- 
tween 6.8 and 9.1 kpc ( Jonker fe Nelemansl 2004) o r be- 



tween 5.9 ± 0.9 and 7.7 ± 0.9 kpc (I Wolff et al.ll2005h . the 



smallest and largest values corresponding to hydrogen rich 
and poor material, respectively. In this paper, we use the 
extreme values of the distance derived so far: 5 and 10 kpc. 

EXO 0748-676 light cur ves exhibit dipping a ctivity and 
eclipses every 3.8 h (e.g. Parmar et al. Il986h . Dips and 
eclipses are due to the central X-ray source being ob- 
scured by some structure above the disk, and occulted 
by the co mpanion star, res pectively, at every orbital pe- 
riod (e.g. iFrank et al.lll987t ). Their presence indicates that 
EXO 0748-676 is vi ewed at an inclination of w 75-83° 
(Par mar et al. I ll986h . i.e. almost from the accretion disk 



plane. 

Here, we report the discovery of triple bursts in 
EXO 0748-676. This is the first time that unambiguously 
and repeatedly triple bursts are detected in an accret- 
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ing neutron star. We inspect EXO 0748-676 non-bursting 
emission to conclude that the accretion rate is likely con- 
stant throughout the analyzed set of XMM-Newton obser- 
vations (Sect. [3]). We derive the properties of the 76 bursts 
detected, using the 5-10 keV light curve to minimize the 
contamination from the dipping activity (Sect. HJ. We 
perform a time-resolved spectral analysis of a sample of 
bursts not affected by dipping. We compare the triple 
bursts with the single and double bursts detected dur- 
ing the same XMM-Newton observations, and also with 
those detected during the EXOSAT observations (Sect [5]). 
We discuss the burst properties in the context of the well- 
studied helium-triggered hydrogen-helium burning ther- 
monuclear flash model and in the context of the poorly- 
studied hydrogen-triggered hydrogen(-helium?) burning 
flash model (Sect[H|). 



2. Observations and data reduction 



The XMM-Newton Observatory (jJansen et al.l l200lh in- 
cludes three 1500 cm 2 X-ray telescopes each with a 
European Photon Imaging Camera (EPIC) at the fo- 
cus. Tw o of the EPIC ima ging spectrometers use MOS 
CCDs dTurner et al. 



20011 ) and one uses PN CCDs 



(IStriider et alj 120011). Reflect ion Grating Spectrometers 



(RGS, Iden Herder et al.ll200lh are located behind two of 
the telescopes. 

EXO 0748-676 was observed by XMM-Newton on sev- 
eral occasions with different instrument configurations. In 
this paper, we use the 7 observations performed between 
September and November 2003, during XMM-Newton 
revolutions 692 to 719, for a total exposure time of 158 
hours (Table Q}. We focus on the EPIC PN data (0.1- 
10 keV) that were all obtained in small window mode with 
the medium optical blocking filter applied. In this mode, 
only a fraction of the central CCD chip, corresponding to 
63x64 pixels or 4.'3x4f4, is read out. This allows a time 
resolution of 5.7 ms to be reached, and significant photon 
pile-up occurs only for count rates ^100 counts s _1 . 

We reduce the observations using versions 5.4.1 to 6.1.0 
of the science analysis software. Only single and double 
events (patterns to 4) are selected. The ratio of back- 
ground to source persistent intensity is generally ^51% 
and reaches occasionally ~10% during some episodes of 
higher background level due to enhanced solar activity. 
We do not discard any data interval for this analysis. For 
the light curves and estimate of rough burst properties 
(peak count rates, number of counts in a burst), we extract 
source events from within a circle of 40" radius centered 
on the position of EXO 0748-676. For the burst spectral 
analysis, we extract source events from an annulus with 
an outer radius of 40" and a 9.3" inner radius, in order 
to avoid pile-up effects which become significant near the 
burst peaks. Such an annulus contains 63% of the events 
from the circular region. 
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Fig. 1. Color (counts in the 5-10 keV band divided by 
counts in the 0.3-5 keV band) of EXO 0748-676 non- 
bursting emission as a function of the soft (0.3-5 keV) 
intensity during the 7 XMM-Newton observations whose 
revolution number is indicated in the right corner. Error 
bars are shown on one third of the points which repre- 
sent 60 s each. The intensity and spectral changes during 
revolution 719 (bursts and eclipses excluded), and hence 
the boomerang track associated to that revolution, can 
be explained by changes in the ionization level and in 
the amount of a bsorbers in front of a constant underly- 
ing X-ray source ( Diaz Trigo et all ioOO). Remarkably, the 
data obtained during the other revolutions (spanning sev- 
eral months) follow the same track. This indicates that 
the spectral changes during these observations can be ex- 
plained by the same phenomenon and that the underlying 
X-ray emission is likely the same during the 7 observations 
studied here. 



3. Non-bursting emission 

The PN light curves of the 7 observations of EXO 0748- 
676 are shown in Fig. lA.ll in the appendix, in the "hard" 
5-10 keV (panel a) and in the "soft" 0.3-5 keV energy 
band (panel b), while panel c shows the color (counts in 
the hard band divided by those in the soft band) as a 
function of time. In addition to X-ray bursts, EXO 0748- 
676 light curves exhibit eclipses every 3.8 h and dipping 
activity. 

Dipping is associated with spectral hardening 
(Fig. IA.il c). the soft light curve (Fig. IA.il b) being sub- 
stantially more affected than the hard one (Fig. IA.1I a) . 
The dipping activity is highly irregular on timescales of 
days, with the soft light curve shape clearly changing from 
one XMM-Newton observation to another. While a dip- 
ping pattern covering a limited phase range of ~0.7-0.9 
is easily recognizable in some observations (e.g. revolution 
693), this is not the case in others (e.g. revolution 692) 
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where the soft intensity and hardness display erratic vari- 
ability at all orbital phases. 

iDi'az Trigo et all (|2006h analyzed the PN data of 
EXO 0748-676 during revolution 719. Excluding eclipses 
and bursts, they define six intensity stages, and call the 
highest range the "persistent" level, and the lowest one 
the "Dip 5" level. They extract one spectrum per level 
and fit them simultaneously with a continuum model af- 
fected by a neutral and by an ionized absorber. While the 
parameters of the continuum are forced to be the same for 
the six spectra, the parameters of the absorbers are left 
indepen dent. By successfully fit ting the spectra with this 
method, iDfaz Trigo et al. (2006) demonstrate that all the 
intensity and spectral changes from EXO 0748-676 during 
revolution 719 can be simply explained by changes in the 
ionization level and in the amount of the absorbers that 
are located inside the system, while the underlying X-ray 
source is staying constant. Not only the "dipping" spec- 
tra are affected by the local absorbers, but also the "per 



sisten t" one (see Fig. 3 and Table 7 of I Diaz Trigo et al 



l2006lV The term "persistent" , if understood as un-affected 
by local absorbers, is therefore un-appropriate in this case, 
EXO 0748-676 being "dipping" (if understood as affected 
by local absorbers) at all intensity levels and at all phases 
during that observation. 

In Fig. [TJ we show the color-intensity diagram of the 
non-bursting emission of EXO 0748-676 during all the 
XMM-Newton revolutions examined here (from 692 to 
719). Remarkably, the data from the different observa- 
tions overlap perfectly in this diagram. Since the inten- 
sity and spectral changes during revolution 719, and hence 
the boomerang-like track associated with that revolution 
in Fig. [TJ could be explained by changing absorbers in 
front of a constant X-ray source, we deduce that the spec- 
tral changes during the other observations, and thus their 
associated boomerang tracks in Fig. [lj can be explained 
similarly and that the underlying X-ray emission is the 
same during revolutions 692 to 719. If the underlying X- 
ray emission had been different from one observation to 
another, the various tracks would have appeared shifted 
with respect to each other, rather than being overlapping, 
unless a change in the underlying X-ray emission had the 
very same spectral signature in the color-intensity diagram 
as the change due to the absorbers, which is unlikely. 

Therefore, in this paper, we assume a constant un- 
derlying X-ray flux, and hence a constant accretion rate, 
for EXO 0748-676 throughout the observations 692 to 719. 
We will consider the unabsorbed X-ray flux derived by 
Diaz Trigo et all (|2006h from revolution 719 as representa- 



tive of the unabsorbed flux for all the observations studied 
here. We note that the maximum PN count rate reached 
by EXO 0748-676 outside bursts varies at most by a fac- 
tor 1.3 from one observation to another, which could be 
considered an upper limit to the variation of the underly- 
ing flux, if the above-mentioned interpretation of changing 
absorbers in front of a constant underlying X-ray emitter 
throughout revolutions 692 to 719 was not correct. 



Us ing the data from revolution 719, iDiaz Trigo et al 
(200(1) derive a persistent 0.6-10 keV absorbed flux of 
2.25 x 10~ 10 erg cm~ 2 s -1 . This corresponds to a 0.6- 
10 keV unabsorbed (without any attenuation) flux of 
2.81 x 10~ 10 erg cm~ 2 s _1 . Introducing a break at 
50 ke V (as found in th e 0.1-100 keV BeppoSAX spec- 
trum, ISidoli et al.l 120051 ) in the power-law component of 
their model, we derive, in the 0.1-100 keV energy band, 
an unabsorbed flux of 8.44 x 10~ 10 erg cm~ 2 s _1 that we 
consider as the representative bolomctric underlying flux 
of EXO 0748-676 for revolutions 692 to 719. At 5 kpc, this 
implies an unabsorbed luminosity of 8.44 x 10 35 erg s _1 
in the 0.6-10 keV band, and 2.52 x 10 36 erg s" 1 in the 
0.1-100 keV band. At 10 kpc, this implies a luminos- 
ity of 3.36 x 10 36 erg s" 1 in the 0.6-10 keV band and 
1.01 x 10 37 erg s _1 in the 0.1-100 keV band. We estimate 
the relative error on the fluxes and luminosities due to the 
spectral fit uncertainties to be ^3%. 



4. Bursts 

76 X-ray bursts are recorded during the 7 XMM-Newton 
observations of EXO 0748-676 (Table UJ). As illustrated in 
Figs.[5]and[n](see also the complete set of bursts in Fig. IA.ll 
of the appendix) , most bursts occur after a wait time since 
a previous burst of typically 3 h, while others occur after 
a wait time of only ~12 min. We will see in Sect. 14.21 
that there is a clear distinction between these long and 
short wait-time bursts. We call several consecutive bursts 
separated by short wait times a bursting event. The obser- 
vations contain events of one, two and three bursts, which 
we refer to as singlets, doublets and triplets, respectively. 
The burst types are denoted by SI for a singlet, Dl and 
D2 for the first and second burst in a doublet, respectively, 
and Tl, T2 and T3 for the first, second and third burst 
in a triplet, respectively. 

The observations contain 52 bursting events: 33 sin- 
glets, 14 doublets and 5 triplets. On average, one event 
occurs every 3 h, one singlet every 5 h, one doublet every 
11 h and one triplet every 20 h. This is the first time that 
triplets are unambiguously and repeatedly detected in an 
accreting neutron star. The bursts from all the triplets are 
shown in Fig. 2] together with the bursts from one singlet 
and from one doublet. 

To reliably compare all the bursts despite the fact that 
many of them occur during dipping (Fig. IA.lj) . unless oth- 
erwise stated, we base the bulk of our analysis on the 
hard 5-10 keV energy band because it is much less con- 
tamina ted by the dipping acti vity than the soft band (cf. 
Sect. 131 IDfaz Trigo et al1l2006l) . In addition, in Sect. [DEI 
we perform a spectral analysis in the full 0.1-10 keV band 
of a sample of bursts that are the least affected by dip- 
ping, from which we derive conversion factors that we use 
to convert the 5-10 keV parameters values into bolometric 
ones in the rest of the analysis. Bursts 26 and 38 contain 
data gaps and are excluded from part of the analysis. 
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Fig. 2. 5-10 keV EPIC PN light curve of EXO 0748-676 with a binning of 60 s during XMM-Newton revolution 693. 
The X-ray bursts are numbered and occur either in singlets, doublets or triplets. 
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Fig. 3. 5-10 keV EPIC PN light curve of EXO 0748-676 
with a binning of 3 s showing the triplet of bursts 13, 14 
and 15. Time is given in minutes from the peak time of 
the burst 13. 

4.1. Parameters definition 

Prior to each burst, we select a reference interval free of 
eclipse, instrumental data gaps or other bursts, in the 1-s 
resolution light curve. We define the burst start time when 
the count rate reaches 3 a above the level in the reference 
segment, as the burst rises, and the burst stop time when 
the intensity drops back to a level of 1.1 times the mean 
intensity in the reference segment, as the burst decays. 
We define the burst duration as the separation between 
these burst stop and start times. The peak time is the 
time at which the count rate is maximum, again in the 1-s 
resolution 5-10 keV light curve. We define the wait time 
of a burst as the separation between its peak time and the 
peak time of the previous burst. We take the event wait 
time to the next event as the separation between the peak 
time of the first burst in the event to the peak time of the 
first burst in the following event. 

We determine the total net counts in a burst by sum- 
ming all the counts detected between the burst start and 



stop times and subtracting from this value an estimate of 
the number of persistent-emission counts that were emit- 
ted over the burst duration (average persistent count rate 
in the reference segment multiplied by the burst duration) . 
The net counts in a burst is a tracer of the burst fluence. 
We define the net counts in the event (tracer of the event 
fluence) as the sum of the net counts of each burst in the 
event. 

We apply a multiplicative factor of 8.1 x 
10 -11 erg cm~ 2 count -1 to convert the 5-10 keV 
net counts into the bolometric fluence, and a factor of 
8.2 x 10~ n erg cm" 2 count -1 to convert the 5-10 keV 
peak count rate into the bolometric peak flux (see details 
in Sect.gXU). 

The a value of a burst is defined as / p x t wa it/Ehi 
where / p is the bolometric persistent flux, t wa i t the burst 
wait time and E\> the bolometric burst fluence. As the 
persistent flux depends on the accretion rate, the prod- 
uct /p x t wait is a tracer of the mass accumulated on the 
neutron star via accretion between the preceding burst 
and the considered one. In the context of thermonuclear 
flash models, the burst fluence is expected to be propor- 
tional to the amount of nuclear fuel available. At accretion 
rates where the fuel does not burn stably outside bursts, 
the burst fluence is thus expected to be proportional to 
the accumulated mass, and a to be constant, in a given 
bursting thermonuclear regime. Here, we determine the a 
values assuming a constant bolometric persistent flux of 
8.44 x 10~ 10 erg cm~ 2 s -1 (Sect. [3]). 

4.2. Parameters distributions 

The peak time, wait time and type of all the bursts are 
listed in Table EU The mean, minimum and maximum 
values of the burst wait time, peak flux, duration, fluence 
and a are given in Table [2] for each burst type and for 
some combinations of types. 



6 



L. Boirin et al.: Burst triplets in EXO 0748-676 



Table 2. Mean, minimum and maximum values of the burst peak flux, duration, fluence and a value, for each 
burst type, derived from the 1 s resolution EPIC PN 5-10 keV light curve, assuming a conversion factor of 8.1 x 
10 -11 erg cm -2 count -1 between the 5-10 keV net number of counts and the bolometric burst fluence, a conversion 
factor of 8.2 x 10 -11 erg cm -2 count -1 between the 5-10 keV peak count rate and the bolometric burst peak flux 
(Sect. I4.6.2[) and a bolometric persistent flux of 8.44 x 10 -10 erg cm -2 s -1 (Sect. [3]). 



Type t W!1 it (h or min) Peak flux Duration (s) Fluence a 

(10 -9 erg cm -2 s -1 ) (10 -7 erg cm -2 ) 





Mean 


Min 


Max 


Mean 


Min 


Max 


Mean 


Min 


Max 


Mean Min 


Max 


Mean 


Min Max 


SI 


3.20 h 


2.23 h 


4.64 h 


7.2 


5.5 


11 


116 


55 


177 


1.9 


0.90 


3.1 


53 


42 


102 


Dl 


2.87 h 


2.04 h 


3.70 h 


7.0 


6.1 


8.6 


100 


65 


156 


1.7 


1.1 


2.2 


52 


45 


61 


Tl 


2.46 h 


2.23 h 


2.79 h 


6.9 


6.4 


7.5 


81 


61 


94 


1.4 


1.2 


1.6 


54 


49 


57 


SI, Dl, Tl 


3.02 h 


2.04 h 


4.64 h 


7.2 


5.5 


11 


109 


55 


177 


1.8 


0.90 


3.1 


53 


42 


102 


D2 


13.20 min 


8.40 min 


19.02 min 


4.7 


1.7 


7.5 


39 


22 


53 


0.61 


0.26 


1.1 


12 


6 


17 


T2 


12.18 min 


9.42 min 


15.36 min 


5.2 


3.2 


7.9 


50 


28 


116 


0.73 


0.38 


1.5 


10 


5 


16 


T3 


11.70 min 


10.62 min 


13.02 min 


2.9 


1.6 


6.3 


23 


15 


42 


0.30 


0.091 


0.87 


35 


8 


63 


D2, T2, T3 


12.66 min 


8.40 min 


19.02 min 


4.5 


1.6 


7.9 


38 


15 


116 


0.57 


0.091 


1.5 


16 


5 


63 


All types 


2.05 h 


8.40 min 


4.64 h 


6.3 


1.6 


11 


86 


15 


177 


1.4 


0.091 


3.1 


40 


5 


102 



Where relevant, we indicate the skewness and kurtosis 
of the samples. The skewness of a sample . . . , Yjv 

is defined as - Y) 3 /(N - l)s 3 , where Y is the 

mean, s the standard deviation, and N the number of 
data points. The skewness is expected to be for a 
symmetric distribution, < for a skewed left one, and 
> for a skewed right one. The kurtosis, defined as 
Eili( y ; " Y) 4 /{N - l)s 4 ] - 3, is expected to be for 
a normal distribution, < for a flatter one, and > for a 
more peaked one. 

We use the two-sample Kolmogorov-Smirnov test to 
estimate the probability, Pg am b c , that two samples a and 
b come from the same distribution. Since there are only 
5 triplets (as compared to 14 doublets and 33 singlets), 
the probabilities involving the triplets are inferred from 
particularly small numbers and should be considered with 
care. 

4.2.1. Wait time 

The histograms of the wait times are shown in Fig. \E\ for 
each burst type and for some combinations of burst types. 
The wait times range between 8.4 min and 4.6 h, but we 
observe no burst with a wait time between 20 min and 
2 h. Two distinct groups of bursts are thus clearly visible: 
one group made up of the SI, Dl and Tl bursts which 
occur after a mean wait time of 3.0 h, and another group 
consisting of the D2, T2 and T3 bursts which occur after 
a mean wait time of only 12.7 min. We will refer to these 
two groups as the bursts with long wait times (LWT) and 
short wait times (SWT), respectively. This distinction is 
confirmed by the very low probability (P sanfc sh ° rt = 7 x 
10 -14 ) that the wait times from the two groups come from 
the same distribution. The SWT histogram and the LWT 
histogram have the same skewness of 0.64 , and a similar 
kurtosis of -0.57 for the LWT and -0.31 for the SWT. 



Inside the LWT group, the wait time distributions 
from two different burst types are not significantly differ- 
ent (P s sl mc D1 = 37%, P S s a 1- C T1 = 23%, P ^ -T1 = 25%). 
Inside the SWT group, the D2, T2 and T3 bursts 
also seem to follow the same wait times distribution 

( P S?~~3 T2 = 96 %, P Ja 2 mo T3 = 7 ^%, P ° 2-T3 = 35%). 

Fig. E] shows the histograms of the event wait 
time to the next event. While the doublets and 
the triplets seem to follow the same distribution 
( p same let_triplct = 58 %), the singlets follow a significantly 
different distribution than that of the doublets and triplets 

( p !kmo Gt " doublct = 0.0015% and P ™mi ct_t " plct = 0- 2 %)- 
For a new event to occur, one has to wait on average 4.2 h 
after a triplet, 3.7 h after a doublet, but only 2.7 h after 
a singlet. 

4.2.2. Peak flux 

The histograms of the burst peak flux are shown in Fig. [7] 
(top; see also Tabled]). Despite their overlap, the LWT and 
SWT peak flux histograms are inconsistent with coming 
from the same distribution (P i™# -short = 1.5 x 10 -8 ). 
The mean peak flux of the LWT bursts is a factor 1.6 
higher than that of the SWT bursts. The mean peak 
flux of the LWT bursts, 7.2 x 10 -9 erg cm -2 s -1 , cor- 
responds to a luminosity of 2.2 x 10 37 erg s -1 at 5 kpc or 
8.6 x 10 37 erg s -1 at 10 kpc. The mean peak flux of the 
SWT bursts, 4.5 x 10 -9 erg cm -2 s -1 , corresponds to a lu- 
minosity of 1.3 x 10 37 erg s -1 at 5 kpc or 5.4 x 10 37 erg s -1 
at 10 kpc. 

The LWT peak flux histogram is very peaked (kur- 
tosis of 5.68) while that of the SWT group is flat (kur- 
tosis of -1.06). The SI, Dl and Tl bursts have similar 
peak flux distributions (P = 84%, P f a 1- c T1 = 96%, 

p ?an7c T1 = 99%). In the SWT g rou P: the D2 and T2 bursts 
have similar peak flux distributions (P ^ l 2 ^ c T2 = 99%), 
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Fig. 4. 5-10 keV EPIC PN light curves of EXO 0748-676 
with a binning of 1 s showing the bursts profiles in one sin- 
glet, in one doublet and in all of the five detected triplets. 
The burst number and its type are indicated in each panel. 
Time is given in seconds from the burst peak time. 



o 



o 
o 
o 

"o 



30 f 



25 



20 



15 



10 



4 - 



D All types 
US1, D1.T1 
': • D2, T2, T3 




2 3 4 5 
Wait time (h) 



10 



15 

Wait time (min) 



20 



S1 

Idi 

! D2 

In - 

] T2 
!T3 



I I 


i i i i i i 

': j D2, T2, T3 Si D2 




T2 ' 








■ t3~ 








I • • 


i : ■ tit 






1 ' ' . 1 








1 *. '-. 



Fig. 5. Histograms of the wait times. The bottom plot 
zooms in on the top one. The histograms are built inde- 
pendently for each burst type and for some combinations 
of types, using always 10 bins. Consequently, the bin size 
is different from one histogram to another. Furthermore, 
the size of the first and last bins is adjusted to match the 
exact observed range of wait times for a given histogram. 



while the T3 bursts seem to follow a marginally differ- 
ent one (P ? a 2 mc T3 = 4.5%, P ^ = 3.6%), with lower 
values. 

In a given doublet (triplet), the peak flux of the D2 
(T2) is significantly lower than that of the Dl (Tl) in 
10/13 (3/5) cases, and similar to it in the other cases (see 
Fig. |4j. The peak flux of the T3 burst is always signifi- 
cantly less than that of the T2 burst, except for the strong 
T3 burst number 24. 

4.2.3. Fluence and duration 

The histograms of the burst fluence are shown in Fig. [7] 
(middle; see also Tabled]). The fluence sample is clearly bi- 
modal (P^e" Bhort = 7 - 6 xl0~ 14 ), the SWT bursts being 
less fluent by a factor ~3 than the LWT bursts on aver- 
age. The mean fluence of the LWT bursts is 1.8 x 10~ 7 erg 
cm~ 2 and corresponds to 5.3 x 10 38 ergs at 5 kpc or 
2.1 x 10 39 ergs at 10 kpc. The mean fluence of the 



SWT bursts is 5.7 x 10~ 8 erg cm -2 and corresponds to 
1.7 x 10 38 ergs at 5 kpc or 6.8 x 10 38 ergs at 10 kpc. 

Inside the LWT group, the flucnces from the SI 
and Dl bursts seem to be drawn from the same dis- 
tribution (Pfa 1 ,™? 1 = 26%), while the Tl bursts would 
follow a marginally different one (P gamtP = 4.1%, 
P samc T1 — 5.4%), the Tl being on average less flu- 
ent than the Dl or the SI bursts. Inside the SWT 
group, the D2 and T2 have similar fluence distributions 
(P = 93%), while the T3 have marginally different 

(lower) ones (P s D a 2 m - T3 = 0.7%, P ™ mc T3 = 3.6%). The his- 
tograms of the burst duration (not shown here) are very 
similar to the fluence ones. 

In a given doublet, the D2 burst is always shorter and 
less fluent than the Dl. In a given triplet, the duration 
and fluence of the T2 is always less or equal to that of the 
Tl. The T3 burst, except for the strong T3 burst number 
24, has a duration and fluence less than that of the T2 
burst. 



L. Boirin et al.: Burst triplets in EXO 0748-676 



Wait time (h) 

3 4 5 




Wait time (h) 

3 4 5 





4 / 

\- 1 -ft-'' 












o Singlet 


w 


Doublet 




X Triplet 



[Persistent flux x wait time] (erg cm ) 



5 10 15 20 

[Persistent flux x wait time] (10 6 erg cm" 2 ) [Persistent flux x wait time] (10 6 erg cm" 2 ) 



Fig. 8. Left: Burst fluence as a function of persistent flux times wait time (accumulated mass tracer). The line indicates 
the mean a value for the LWT bursts (a — 53). Middle: zoom in on the left panel. Right: Event flucncc as a function 
of the accumulated mass tracer. The line is the same as in the left panel. 
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Fig. 6. Histograms of the event wait times to the next 
event, built independently for the singlets, doublets, 
triplets and for all events, using 8 bins. 



4.2.4. Alpha 

The histograms of the a values are shown in Fig. [7] 
(see also Table The sample is clearly bimodal 
( p iame" sh ° rt = 1-0 xlCT 12 ), the mean a value being 
smaller by a factor ^3.3 for the SWT bursts (a = 16) 
than for the LWT ones (a — 53). The a histogram for the 
SI bursts is very peaked (kurtosis of 11.4) compared to all 
the other burst types (kurtosis between -0.4 and -2). Inside 
the LWT group, the SI, Dl and Tl bursts seem to follow 
the same a distribution (P = 90%, P j^P = 48%, 

PPamJ 1 = 60%). Inside the SWT group, the T3 bursts 
seem to follow a marginally different a distribution from 



the D2 and T2 bursts (P 
while P? a 2 m - T2 = 32%). 



D2-T3 
same 



1.9%, P ""J™ = 3.6% 



4.3. Relation between burst fluence and accumulated 
mass 

Fig. [8] (left and middle panels) shows the burst fluence as 
a function of the product of the persistent flux times the 
wait time, which is a tracer of the mass accumulated by 
accretion before the burst. Here, the accumulated mass 
tracer is simply proportional to the wait time because we 
assume the underlying persistent flux is constant (Sect. [21) ■ 

Among the LWT types, the fluence is strongly corre- 
lated with the accumulated mass tracer, with correlation 
coefficients of 0.91 for the SI, 0.89 for the Dl and 0.81 
for the Tl bursts. For each LWT type, the fluence is well 
described by a linear function of the accumulated mass 
tracer (E^ — or 1 x [/ p x i wa i t ] + Eq). Eq is the offset en- 
ergy at zero burst interval. The best-fit parameter values 
a and Eq are listed in Table [3] Within the 90% confidence 
range, both a and Eq are consistent with being the same 
for the SI, Dl and Tl bursts. For all the LWT bursts to- 
gether, the a value obtained from the fit is 46 ± 4 and 
E is (—2.3 ± 1.7) x 10~ 8 erg cm~ 2 . Eq is consistent with 
within 3rr, and consistent with being negative within 
the 90% confidence range (1.65er). A negative value of Eq 
would indicate that a part \Eq \ — 2.3 x 10~ 8 erg cm -2 of 
the total available energy is not burned in the burst. This 
corresponds to an energy of 6.9 x 10 37 ergs at 5 kpc or 
2.7 x 10 38 ergs at 10 kpc, equivalent to 13% of the mean 
energy released by the LWT bursts. 

Among the SWT types, only the T3 bursts show a 
strong correlation (coefficient of 0.86) between the flu- 
ence and the accumulated mass, which can be described 
by a line with best-fit parameters a = 1.2 ± 0.4 and 
Eq — (— 4.5±4.6)xl0 -7 ergcm~ 2 (90% confidence errors). 
The flucncc of the D2 or T2 bursts is much less strongly 
correlated with the accumulated mass, the coefficients be- 
ing 0.53 for the D2 and 0.62 for the T2, and no accept- 
able linear fit could be obtained. So, apart maybe for the 
T3 bursts, the fluence of a SWT burst does not strongly 
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(middle) and a values (bottom) built as in Fig. [5] 



depend on the amount of material freshly accreted just 
before the SWT burst. 

Combining SWT and LWT types (e.g. the whole sam- 
ple of bursts, or the Dl and D2 combined), we could not 
obtain acceptable (with reduced x 2 <2) linear fits be- 
tween the fluence and the accumulated mass tracer. Thus, 
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Fig. 9. Bolomctric burst fluence as a function of the bolo- 
metric peak flux. The right and top axis indicate the 
corresponding energy and peak luminosity assuming an 
isotropic emission at a distance of 10 kpc. 



the SWT bursts are inconsistent with following the same 
fluence- accumulated mass relation as the LWT bursts. 

Fig. [S] (right) shows the event fluence as a function 
of the accumulated mass tracer. At a given accumulated 
mass, a triplet is generally more fluent than a doublet 
which is more fluent than a singlet. When all the events 
are considered together, the event fluence is only weakly 
correlated to the accumulated mass tracer (correlation co- 
efficient of 0.38). The fluence of the doublets correlates 
with the accumulated mass. However, this correlation is 
not better than the one observed for the Dl bursts only. 
Similarly the correlation is not improved by considering 
the triplets (coefficient of 0.44) rather than the Tl bursts 
(coefficient of 0.81). This effect is visible when the left 
and right panels of Fig. [S] are compared: the dispersion 
is higher for the doublets than for the Dl and higher for 
the triplets than for the Tl. This suggests that the flu- 
ence of a D2 burst (and analogously of a T2-T3 pair) does 
not strongly depend on the amount of material accreted 
before the Dl (or Tl) burst of the considered event. 



Table 3. Best- fit parameters a and Eq, and their 90% 
confidence errors, obtained from fitting E^, as aT 1 x [/ p x 
iwait] + E Q (see text). 



Type 


a 


E (10" s erg cm' 2 ) 


SI 


45 ±5 


-2.5 ±2.3 


Dl 


49 ±11 


-1.1 ±3.7 


Tl 


50 ±40 


1 ± 13 


SI, Dl, Tl 


46 ±4 


-2.3 ± 1.7 
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Fig. 10. Average profiles for each burst type, obtained re- 
spectively from 32 SI, 13 Dl, 14 D2, 5 Tl, 4 T2 and 4 T3 
bursts. The strong T2 and T3 bursts number 35 and 24 
are excluded. The binning is 3 s. Note that the logarith- 
mic scale of the intensity axis makes an exponential decay 
appear as a straight line. 



4.4. Relation between fluence and flux 

Fig. O shows the burst fluence as a function of the burst 
peak flux. While the two quantities are not strongly corre- 
lated for the SI, Dl and Tl bursts (with linear correlation 
coefficients of 0.04, 0.29, and -0.27, respectively), the flu- 
ence and flux are strongly correlated in the case of the D2, 
T2 and T3 bursts, with correlation coefficients of 0.93, 0.94 
and 0.99, respectively. 

4.5. Profiles 

Fig. Uni shows the profiles obtained by averaging the 5- 
10 keV light curves of all the bursts of a given type, except 
the exceptionally strong T2 and T3 bursts number 35 and 
24. After subtracting the persistent count rate level before 
the bursts, the decay of all these profiles is well described 
by a model consisting of two exponential parts. The best- 
fit values of the characteristic decay times, T\ and T2, and 
of the transition time t\ are given in Tabled] The first part 
of the decay is always faster than the second (t\<t-i). The 
decays of the average SI, Dl and Tl profiles are similar 
to each other, with t\ ~ 15 s and t-i ~ 52 s. The D2 
and T2 average profiles decay more rapidly, with both t\ 
and T2 being shorter. For the T3 average profile, only one 
exponential component is detected (including a second one 
does not improve the fit quality) . Its decay time is similar 
to that of the first component of the D2 and T2 average 
profile. 

The top panel of Fig. [TT] shows the average SI, Dl, Tl 
and D2 profiles obtained from bursts having a similar peak 
intensity in the range 77-96 counts s , and a similar wait 
time in the range 2-2.8 h for the SI, Dl and Tl. The best- 
fit model of their decay with two exponential parts is given 
in Table O While the SI, Dl and Tl profiles are similar, 
the D2 profile, at the same peak intensity level, shows a 
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Fig. 11. Top: Average SI, Dl, Tl and D2 profiles, ob- 
tained respectively from 5 SI, 5 Dl, 3 Tl and 3 D2 bursts 
having a similar peak intensity (77-96 counts s _1 ), and a 
similar wait time for the SI, Dl and Tl (2-2.8 h). Time 
is given from the burst peak time and the binning is 3 s. 
Bottom: The strong T3 number 24 compared with the 
SI and D2 average profiles from the top panel. 

different and much faster decay which can be modeled by 
only one exponential decay. 

The only T2 and T3 bursts that have a peak intensity 
in the range 77-96 counts s _1 are the exceptionally strong 
ones number 35 and 24, respectively. The best-fit model 
of their decay is given in Table The decay of the strong 
T2 burst number 35 shows a fast part followed by a slow 
flat component. Interestingly, the decay of the strong T3 
burst (bottom panel of Fig. [TTj) is very similar to that of 
the average D2 profile at the same peak intensity level. 

4.6. Spectral properties 
4.6.1. Color intensity diagram 

Fig.[T2lshows the color, C, as a function of the hard inten- 
sity, /, during the decay of the averaged profiles obtained, 
respectively, from 15 SI, 8 Dl, 10 D2, 2 Tl, 1 T2 and 
1 T3 bursts chosen for being the least contaminated by 
the dipping activity. We find that, for each burst type, 
log(C) can be described by a linear function of log(J) 
(log(C) = a x log(7) + b). For a burst emitting like a 
blackbody, the bolometric flux F is expected to depend 
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Table 4. Best-fit parameter values and 90% confidence errors for the decay of the average profiles by burst type 
(Fig. [10]). We subtract the persistent count rate level before the burst and fit each persistent-subtracted burst decay 
using a continuous function of two exponential parts: Boe~ l l Tl between to and ij., and i? 1 e~ t / 7 ' 2 between t\ and ti. 
We set to to the burst peak time (0 s), t% to the first time where the persistent-subtracted intensity, B, reaches a level 
of 0, and B\ to S e' _ * 1 ^' ri - ) . The free parameters are the persistent-subtracted intensity, B , at the burst peak time, 
the characteristic decay times, n and T2, and the transition time t±. We indicate the time, i s top, where the intensity, 
/, reaches back a level of 1.1 times the level before the burst. For the profiles whose fit was not improved by including 
a second exponential component, we give only the decay time of the first and only exponential component. 



Type 


B 

(counts s~ ) 


Tl 

00 


T'2 

00 


ti 
00 


^stop 

00 


SI 
Dl 
Tl 
D2 
T2 
T3 


' a -' -0.5 

74 1+3-0 

71 1+ 4S 
' '-•'--is 

53.9+f'g 
50.51*;? 

i 7 q+8.2 
1 ' "J-0.4 


i5.9±S:i 

13-8l ;°3 

n-olSI 

Q 0+ 11 

10.8l 8 'i 


56.9lgi 
51+ 4 
48+ 2 2 
19+ 2 3 
14+i° 


20.7 +0 2 7 
21.6+^ 
20.6l°; 6 3 
25+1 
15±3 


180 
126 
103 

52 
45 
33 



Table 5. Same as Table for the average profiles by type from bursts having a peak intensity in the range 77- 
96 counts s" 1 and a wait time in the range 2-2.8 h for the SI, Dl and Tl (Fig. fTTjl . 



Type Bo n r 2 ti t stop 

(counts s _1 ) (s) (s) (s) (s) 

SI 70l? 14.0l^ 49l^° 19.9±i;| 116 

Dl 8llf 14.2l ;l 40 ±i 20.3l^ 105 

Tl 68l^ 13.7lJ;l 4 5l2 4 21.01!^ 101 

D2 7511° 11.6± ;| 50 

T2 78lf 13.0±g;| 149<- g 277 a 30.7±^;g 108 

T3 92±g IQjlpj — — 37 



a The second component of this decay is effectively flat, 

which shows in our fits as a long positive or negative exponential decay time. 



only on the area, A, of the emitting region and on the 
blackbody temperature, T, as F cx AT 4 . Since C is a 
tracer of T and / is a tracer of F, any difference in the 
slope a between burst types would indicate a difference in 
the flux to temperature dependence (F oc T 4 ), while any 
difference in the intercept b would indicate a difference in 
the area A of the emitting region. We find that, within 3<r, 
the slope a is consistent with being the same for all burst 
types. So the comparison of the averaged color-intensity 
diagrams by burst types does not indicate differences by 
type in the burst temperature to flux dependence. The 
best-fit values of the intercept b are consistent with being 
the same within 3(7 for the SI, Dl, Tl, T2 and T3 types 
on one hand, and for SI, Dl, D2, T2 and T3 types on the 
other hand, indicating a difference only between the Tl 
and D2 types. This difference suggests that the apparent 
emitting area (or blackbody radius) would be on average 
larger for the Tl than for the D2 bursts. 

Fig. H3] shows the color at the burst peak as a func- 
tion of the intensity at the peak, for the individual bursts 
the least affected by dipping activity. The two quantities 
are strongly correlated. This implies that the peak tem- 
perature distributions have similar properties as the peak 



fluxes distributions shown in Fig. [71 with a clear distinc- 
tion between the SWT and the LWT bursts, the latter 
having higher and less variable peak temperatures. In ad- 
dition, the correlation in Fig. [T3] is consistent with that 
obtained from the averaged decays (Fig. [T2")) . So the re- 
lation between the temperature and the flux is consistent 
with being the same for all burst types and for any time 
of the burst decay including at the peak. 

This analysis shows that there is no indication for 
strong spectral differences between the various burst 
types, the flux or apparent emitting area being approx- 
imately the same at a given temperature during the de- 
cay. There is however marginal evidence that the apparent 
emitting area is on average larger for the first bursts in 
triplets than for the second bursts in doublets. This could 
point to differences in the structure or composition of the 
neutron star atmosphere/outer lay ers (see Sect. [5] and the 
discussion in lGottwald et al.lll987t ). 

4.6.2. Time-resolved spectral modeling 

We select five singlets (bursts 1, 12, 16, 33 and 37, see 
Appendix), one doublet (bursts 4 and 5) and two triplets 
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Fig. 12. Color as a function of hard intensity during the 
decay of the averaged profiles by type obtained, respec- 
tively, from 15 SI, 8 Dl, 10 D2, 2 Tl, 1 T2 and 1 T3 
bursts which are the least affected by the dipping activ- 
ity. Error bars are shown on one third of the points which 
represent each 1 s. 



(bursts 22, 23, 24 and 34, 35, 36) that are the least contam- 
inated by dipping activity and perform their time-resolved 
spectral analysis using version 11.3.1 of the XSPEC pack- 
age. We divide each burst into 3 to 6 intervals, with the in- 
terval around the burst peak chosen shorter than intervals 
in the tail. We extract a spectrum in each burst interval 
from which we subtract a reference persistent spectrum 
obtained from segments just before and after the burst. 
We add a 2% systematic error to account for calibration 
uncertainties. 

Since the net emission of Type-I X-ray b ursts is gener 



ally well described by a blackbody (see e.g. ISwank et al 



1977t iLewin et all 1 19931 ). we model the net 0.1-10 keV 
spectrum with an absorbed blackbody. We fix the hy- 
drogen column d ensity, TVh, of the neutral absorber to 
0.1 x 10 22 cm- 2 (ISidoli et al.l l2005h . We also include a 



Gaussian with a centroid energy set to to describe emis- 
sion below ~2 keV which is in excess of the blackbody 
model. This soft excess may be related to a temporary in- 
crease, due to the burst irradiation, of the ionization level 
of the local absorbing material in the line-of-sight. The 
material becomes more transparent, especially at low en- 
ergy, than it was just before the burst. Consequently, the 
subtraction of the spectrum extracted prior to the burst 
from the spectrum extracted during the burst yields to an 
apparent soft excess in the net spectrum. This component 
will not be investigated further here. A soft excess was re- 
ported in other bur sts from EXO 0748-676 and interpreted 
in a similar way by lAsai fc Dotanl ( 2006^ 1 . 

We obtain good fits for each interval. Fig. Q3] (left) 
shows the derived blackbody temperatures and bolomet- 
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Fig. 13. Color as a function of the hard intensity at the 
burst peak for the individual bursts the least affected by 
dipping activity. 



ric fluxes as a function of time. Fig. [TH (right) shows the 
blackbody radius (top) and the bolometric flux (bottom) 
as a function of the blackbody temperature. For a given 
burst, we average the best- fit parameter values over the 
intervals near the peak on one hand, and over the intervals 
in the tail on the other hand. In Table [HJ we give the range 
of these peak and tail values covered by our burst sample. 
The blackbody surface temperature fcTbb lies in the inter- 
val (0.93 — 2.13) ±0.11 keV, where higher values are found 
in the burst peak and lower values in the burst tail. There 
is no evidence for photospheric expansion. The normal- 
ization of the blackbody component and hence the radius 
of the emitting region is consistent with being constant 
throughout a burst and consistent with being the same 
for all the bursts studied. The average radius of the black- 
body component is 2.4 km assuming a source distance of 
5 kpc, and 4.9 km assuming a distance of 10 kpc, with 
an average uncertainty of 44%. This blackbody radius, 
since derived from a measured color temperature rather 
than the effective temperature, is an undere stimate of the 
physi cal radius of the neutron star (e.g. lLondon et al. 
11984) . 

Using the previous spectral models, we further cal- 
culate the unabsorbed bolometric flux of the blackbody 
component for each interval, and finally the bolometric 
burst fluence by integrating over the whole burst. We 
find that, for our sample of bursts, the bolometric flu- 
ence linearly correlates with the 5-10 keV net counts in 
the burst. This means that the ratio between the bolo- 
metric fluence and the 5-10 keV net counts in the burst 
is consistent with being the same for all bursts types. We 
derive a value of (8.1 ± 0.3) x 10~ n erg cm~ 2 count" 1 for 
this ratio. Analogously we find a ratio of (8.2 ± 0.3) x 
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Fig. 14. Results of the time-resolved spectral analysis of bursts number 1, 12, 16, 33 and 37 (SI) , 4 (Dl), 5 (D2), 22 
and 34 (Tl), 23 (T2), and 36 (T3). Left: Blackbody temperature (top) and bolometric flux (bottom) as a function of 
the time since burst onset. Right: Blackbody radius (top) and bolometric flux (bottom) as a function of the blackbody 
temperature. 



peak flux and the 5-10 keV peak count rate. These ra- 
tios serve to estimate the bolometric fluence and peak flux 
from the 5-10 keV counts for all other bursts, even when 
they are contaminated by dipping. 

5. Comparison with the EXOSAT observations 

In 1985, EXOSAT observed EXO 0748-676 at different 
persistent 0.1-20 keV fluxes be t ween 3 and 18 x lO~ 10 erg 



cm-- s _i (|Parmar et al.lll986t iGottwald et al]|l986t ). 26 



temperature increased, and equivalently (iv) their appar- 
ent blackbody r adius at a given temperature increased. 
Gottwald et al. (1986) hypothesized that the persistent 
flux dependent variations in the burst properties could 
be caused by the flashes changing from helium-dominated 
at high accretion rates, to hydrogen-triggered hydrogen- 
helium flashes at low accretion rates. They further specu- 
lated that double bursts could be a feature of hydrogen- 
triggered hydrogen- helium flashes. 



bursts were observed, including four doubl ets with a burst 
separation of the order of 10 20 min (Gottwald et al. 
1986). The doublet phenomenon only occurred when the 



persistent 0.1-20 keV flux was below 5 xlO" 10 erg cm" 2 ([Gottwald et al.l 119871 ) 
s _1 , and was part of an overall correlation of the burst 
properties with changing mass accretion rate. As the per- 
sistent flux increased, (i) the wait time of the LWT bursts 
increased from ~1.8 up to 16 h, (ii) their shape changed 
from a "slow" (long tail) to a "fast" profile sometimes 
showing photospheric expansion, (iii) their flux at a given 



In 1986, EXOSAT caught EXO 0748-676 at a per- 
sistent 0.1-20 keV fl ux of 5 xlO" 10 erg cm" 2 

11 bursts were observed 



s 

and 

showed a regular pattern with a long recurrence time al- 
ways followed by a short one, reminiscent of the double 
burst phenomenon, but with longer wait times, in the 
range 20-70 min for the D2-like bursts. All the burst prop- 
erties were otherwise consistent with that of bursts in the 
slow mode (i.e. low states) of 1985. 
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Using the spectral model bv iDfaz Trigo et aT ( 20061 ) 
(see Sect. [3]), we derive an unabsorbed flux of 4.7 x 
1CT 10 erg cm" 2 s" 1 in the 0.1-20 keV band for EXO 0748- 
676 during the 2003 XMM-Newton observations, indicat- 
ing that the LMXB was again in a low state. We further 
suggest that the source state in 2003 is likely closer to the 
lowest states of 1985 than to the low state of 1986 since 
the unabsorbed EXOSAT fluxes are probably underesti- 
mated compared to the XMM-Newton ones because of the 
local neutral and ionized absorbers included in the XMM- 
Newton spectral models. 

The comparison of Fig. [TJ] with Fi gs. 6 and 7 of 
Gottwald et al. ( 19861 ) and Figs. 4 and 5 of lGottwald et al 



(|l987l) ndicatcs that the spectral properties of the XMM- 



Newton bursts are all consistent with that of EXOSAT 
bursts in the slow mode (low states). In particular, the 
XMM-Newton values of the blackbody radius (4.9±2.1 km 
on averaged) agree well with the EXOSAT values ob- 
tained in the low states whereas larger radii of ~9 km 
had been obtained in higher states. These variations in 
the blackbody radius for a given temperature are more 
likely related to a change in the structure and composi- 
tion of the neutron star outer layers or atmosphere which 
distorts the blackbod y spectrum than to real variations 
of the emitting area (jGottwald et al.lll987l ). In any case, 
the blackbody radius, since derived from a measured color 
temperature rather than the effective temperature, is an 
under estimate of the phy sical radius of the neutron star 



(e. g. London et al. 1984h . The latter one was estimated 
bv lOzell (|2006l ) to be 13.8 ± 1.8 km in EXO 0748-676 



Fig. US] (top) shows the EXOSAT burst wait time dis- 
tributions together with the XMM-Newton one. While 
no burst at all was detected with a wait time between 
~20 min and ~2 h neither in 1985 nor in 2003, this gap 
was partly filled in 1986 by the D2-like bursts with 20- 
70 min wait times. Fig. [T5] (bottom) displays the EXOSAT 
wait time distributions as a function of the persistent 0.1- 
20 keV flux. There is a global shift of the histograms to- 
wards higher wait times as the persistent flux increases. 
The XMM-Newton wait time distribution, and especially 
the presence of doublets with ^12 min burst separa- 
tion, is consistent with that obtained by EXOSAT when 
EXO 0748-676 was in its lowest state in 1985. 

These comparisons confirm the above-mentioned de- 
pendencies of the LWT bursts properties (shape, spec- 
tral characteristics, wait times) on the persistent flux. 
The presence of SWT bursts could also depend on the 
persistent flux in the following scheme. At low persis- 
tent flux (in the low states of 1985 and in 2003), dou- 
blets are emitted with a burst separation of ^10-20 min. 
Triplets occur in that regime as well. They could have been 
missed by EXOSAT by chance since the satellite looked 
at EXO 0748-676 in the low state only 3 times -9 h, while 
XMM-Newton detected a triplet only once every 20 h. In 
1986, the D2-like bursts have longer wait times of ~20-70 
min. This could be related to a slightly higher persistent 
flux then, although the derived fluxes in the low states of 
1985, 1986 and 2003 are too close to each other to con- 



Table 6. Averaged results of fits of the 0.1-10 keV spectra 
extracted during intervals of 13 bursts, with a model con- 
sisting of an absorbed blackbody and a Gaussian. The pa- 
rameters are the hydrogen column density iVu, the black- 
body temperature fcTbb and normalization Abb (defined 
as -Rj^m/^io where i?km is the source radius in km and dio 
the distance to the source in units of 10 kpc), the Gaussian 
centroid energy E, width a and normalization Nq. For a 
given burst, we averaged the best-fit values from the inter- 
vals near the peak on one hand, and from the intervals in 
the tail on the other hand. A range indicates the spread 
in the obtained averaged values for the different bursts. 
Uncertainties are the average errors of the parameters at 
la confidence level. Xreduccd 1S calculated using 50 degrees 
of freedom. Its average and root mean squared are given. 



Peak 



Tail 



iV H (10 22 cm" : 
kT hh (keV) 
iVbb 

E (keV) 
a (keV) 
JVg (cts cm _2 s 

Xrcduccd 



0.1 (fixed) 0.1 (fixed) 

(1.54 - 2.13) ± 0.12 (0.93 - 1.39) ± 0.09 

(6 - 49) ± 5 (7 - 39) ± 5 

(fixed) (fixed) 

(0.5-0.8) ±0.6 (0.40 -0.84) ±0.10 

(0.15 - 1.94) ± 0.10 (0.08 - 0.36) ± 0.05 

1.1 ±0.2 1.2 ±0.3 



sider this link as certain. At much larger fluxes (in the 
intermediate and high states of 1985), there are no bursts 
with short wait times anymore. 

6. Discussion 

Our findings may be summarized as follows. XMM- 
Newton observed EXO 0748-676 in a low state 7 times, for 
an average of 22 h per observation, between September 
and November 2003. 76 sub-Eddington bursts were de- 
tected, either in single, double or triple events, respectively 
on 33, 14 and 5 occasions. The separation between two 
events is ~ 3 h while the separation between two bursts 
in a doublet or a triplet is only ^12 min. The primary 
bursts of different events are similar in terms of duration, 
peak flux and fluence. Their fluence is strongly correlated 
with the amount of material accreted before the burst. 
The profiles of these bursts show an initial fast decay and 
a second slower decaying component. Quite different are 
the secondary and tertiary bursts. They have a shorter 
duration, lower peak flux and lower fluence. Their flu- 
ence is not strongly correlated with the (little) amount 
of material accreted before the burst, except maybe for 
the tertiaries. The profiles of the secondary and tertiary 
bursts show only the fast decay and lack the second, more 
slowly decaying, component seen in the primaries. There 
is no strong spectral differences between the various burst 
types. At a given accumulated mass, the total fluence of a 
triplet is generally larger than that of a doublet which is 
larger than that of a singlet. For a new event to occur, one 
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Fig. 15. Top: EXO 0748-676 burst wait time distributions 
re corded in 1985 and 1 986 by EXOSAT (from Table 2 
of IGottwald et all Il986l and Table 1 of iGottwald et"aH 
19871) and 2003 by XMM-Newton (this paper). Bottom: 
EXOSAT burst wait time distributions as a function of 
the persistent 0.1-20 keV flux of E XO 0748-676 indicate d 



in the plot (from Tables 1 and 2 of IGottwald et al.lll986l ) 



Note the broken wait time axis in both panels. 



has to wait longer after a doublet or a triplet than after a 
singlet. 

6.1. Comparison with theory: LWT bursts 

How do these results fit within the current theoretical 
framework? Let us first consider the recurrence times and 
energetics of the LWT bursts. There are two ways to make 
bursts (or burst "events" ) with a recurrence time as short 
as 3 hours: mixed H/He ignition at an accretion rate close 
to 10% Eddington, or unstable hydrogen ignition at an 
accretion rate of ~ 1% Eddington. 



6.1.1. Mixed H/He ignition 

The best-studied of these is the first, unstable helium ig- 
nition in a hydrogen-rich environment. At high accretion 
rates, the accumulating fuel layer is hot enough that the 
hydrogen burns via the beta-limited a nd thermally stable 
hot CNO cycle (|Hovle fe Fowlerlll965fh When the helium 



ignites, it burns in a hydrogen rich environment, prolong- 
ing the nuclear energy release and leads to long burst tails. 

The shape of the burst lightcurves from EXO 0748- 
676 is similar to those observed from GS 1826-24 which 



Galloway et al.1 (|2004h argued is bursting in this regime. 



The observed a value supports this interpretation, since 



the typical energy release Q 



4 MeV per nu- 



cle on in the rp-process ( e.g., Schatz et al. 19991 : Table 10 
of Wooslev et all 120041) gives an expected a (e.g., 
Galloway et al.ll2004h Y close to the observed value a — 53: 



2 

c z 



= 48 



4 MeV\ 

Qnuc / 



1.5 



(1) 

where z is the gravitational redshift, and £f, )P are fac- 
tors to account for anisotropy in the burst and persis- 
tent luminosities. They are defined by Lb tP = AircP^b.pfb^, 
where L^p and fb, p are the intrinsic luminosities and 
the observed fluxes, and d the distance to the source 
(|Lapidus fe Sunvaevlll985l : lFuiimotolll988l f1. £& >p > 1 in- 
dicates that the observed flux is lower than it would 
be in the absence of anisotropy. Here we use the esti- 
mates ^r 1 = 0.5+ Icosi and Cr 1 = 2 I cos i I , obtained 

i 1 j 1\ ii/ 

bv lFuiimotol (1988) for a system viewed under inclination 

i whose geometry enhances the emission in the direction 
perpendicular to the disk plane, because of scattering in 
the disk, more strongly when the incident photons come 
from the inner disk (during persistent emission) than when 
they come from the neutron star sur face (during bursts) . 
Using i « 75-83° for EXO 0748-676 (IParmar et alJll986h 
gives £b ~ 1.3-1.6, £ p « 1.9-4.2, and a ratio £ p /£b ~ 1.5- 
2.6. We note that the £& iP estimates could significantly dif- 
fer if other system geometries (e.g. including disk warping) 
were assumed. 

In multizone X-ray burst simulatio ns of mixed H/He 



ignition bursts, I Wooslev et al.1 |200J) found recurrence 
times close to 3 hours for accretion of solar metallic- 
ity material at M = 1.75 x 10~ 9 Mq yr _1 , or lumi- 
nosity M(GM/R) » 2 x 10 37 erg s" 1 for a 1.4 10 



km n eutron star (models zM and ZM of IWooslev et al 
l2004h . This agrees quite well with the observed X-ray 



luminosity from EXO 0748-676, once the anisotropy pa- 
rameter is included. The unabsorbed 0.1-100 keV flux 
is 8.4 x 10 -10 erg cm~ 2 s _1 , giving a luminosity Lx = 
0.48 x 10 37 erg s" 1 (£ p /1.9) for a distance of 5 kpc and 
Lx = 1.9 x 10 37 erg s _1 (£ p /1.9) for a distance of 10 kpc, 
where we normalize the anisotropy parameter to i = 75°. 
Therefore we find that the recurrence times and a val- 
ues of the LWT bursts agree well with the models of 



1 We note that anisotrop y of t h e bur st luminosity was not 
included in the analysis by I6zell l|2006h . The EXO 0748-676 
neutron star mass and radius determinations in that paper 
are likely unaffected by the anisotropy factor, as long as it is 
constant through the burst, since they come from taking the 
ratio of the peak flux and the flux in the tail. However, the 
lower limit on the source distance, d > 9.2 ± 1.0 kpc should be 
reduced by a factor ^ 2 , or 15-30%, giving d > 7.3- 8.1 kpc 
for £(, = 1.3-1.6, using the estimate of lFuiimotol l|l988h . 
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Wooslev et al.l (|2004h . particularly when the anisotropy 
factor is included. 

However, the total energies of the LWT bursts in 
EXO 0748-676 are somewhat lower than in the mod- 
els, or than observed for GS 1826-24. In the mixed 
H/He regime , the ignition column dept h is j/i gn ~ 2 x 
10 g cm" " 2 (jCumming fc Bildstenl 12000) giving a burst 
energy E hurst = 47ri? 2 Q nuc y ign /(l + z) = 7.4 x 10 39 ergs 
(assuming complete burning of the whole surface, ra- 
dius R = 10 km, and energy release Q n uc = 4 MeV). 
This estimate agrees well with models zM and ZM of 



Wooslev et al.l ([2004D (see their Table 9) and with the 



observed burst energies in GS 1826-24 (5.3 x 10 ergs 



at 6 kpc, iGallowav et al.l 120041 ) . while the mean burst 
energy for the LWT bursts in EXO 0748-676 is (0.5- 
2.1) x 10 39 ergs for the distance range 5-10 kpc. Note that 
anisotropy, which wou ld tend to decrease th e burst lumi- 
nosity for GS 1826-24 (Gallow ay et al.ll2004l) and increase 
it for EXO 0748-676 (jFuiimotol Il988l ). is not included in 
these estimates. 

The same is true for the peak luminosity. The high 
base temperatures > 10 9 K reached in mixed H/Hc 
bursts give a peak luminosity > 10 38 erg s" 1 (e.g., 
Wooslev et al.ll200"4h , approaching the Eddington luminos- 
ity. Ignoring anisotropy, in GS 1826-24, the bursts peak 
flux was 3 x 10" 8 erg cm" 2 s" 1 , giving L p ea k = 1.3 X 
10 38 erg s -1 at 6 kpc (jGallowav et all 120041) . In contrast 



the mean peak flux for the LWT bursts in EXO 0748- 
676 is 7.2xl0~ 9 ergs cm -2 s -1 , giving L pca k = (0.22- 
0.86) x 10 38 erg s" 1 for the distance range 5-10 kpc. 
Independent of the assumed distance, it is a factor of 7 
below the p eak flux of t h e brig ht radius expansion burst 
observed bv IWolff et all (l2005h from EXO 0748-676, fur- 
ther showing that these LWT bursts are faint relative to 
the Eddington luminosity. 



6.1.2. Hydrogen ignition 

The low burst energies and peak luminosities suggest a 
different explanation for the LWT bursts from EXO 0748- 
676, that they are triggered by unstable hydrogen ignition. 
Hydrogen burning is unstable for temperatures less than 
« 8 x 10 7 K, when the CNO is no l onger beta-limited, an d 
can trigger thermonuclear flashes (jFuiimoto et al. I ll98lh . 

Bildsten ( 19981) estimates that unstable hydrogen ig- 
nition occurs for accretion rates M < 2 x 10" 10 M yr" 1 , 
or luminosities Lx < 2 x 10 36 erg s" 1 . This matches the 
luminosity of EXO 0748-676 if the source is located at the 
closer end of its distance range, and the anisotropy fac- 
tor is small. The observed recurrence time of « 3 hours 
matches the recurrence time expected close to the transi- 
tion between unstable and stable hydrogen burning. The 
maximum temperature at which hydrogen can unstably 
ignite, close to 8 x 10 7 K, corresponds to a column depth o f 
approximately 10 7 g cm" 2 (see Fig. 1 of lCummine 2004 ). 
At an accretion rate of 1% of the Eddington rate, this 
column is accreted in 3.2 hours. 



There ar e two possible outcomes of unstable hydro- 
gen ignition (|Fujimoto et al ll98ll IPeng et al.ll2006f ) . If the 
hydrogen ignition depth is > 5 x 10 7 g cm -2 , the mini- 
mu m column dep th at which helium can ignite (Fig. 1 
of ICumming 2004 ). the increase in temperature following 



hydrogen ignition triggers ignition of helium by the triple 
alpha reaction, and a mixed hydrogen/helium flash occurs. 
At smaller ignition depths, which occur near the transition 
between unstable and stable hydrogen burning, the hydro- 
gen flash is not able to trigg er helium ignition. This case 
has recently been modeled bv lPeng et al. ( 20061 ). They in- 
clude, for the first time, sedimentation of heavy elements 
in the accumulating fuel layer. For an accretion rate close 
to 1% of the Eddington rate, they find that the hydrogen 
flash reaches a peak luminosity of 5 times the accretion lu- 
minosity, or w 0.1 x 10 38 erg s — 1 . This is lower than the ob- 
served peak luminosities of the LWT bursts by a factor of 
a few. Burning a column depth of 10 7 g cm -2 of hydrogen 
to helium gives an expected burst energy 0.6 x 10 39 ergs 
(using the energy release Q n uc = 6.0 x 10 18 erg g" 1 appro- 
priate for the hot CNO cycle). This is within the range of 
observed energy for the LWT bursts. 

Although the peak luminosity found by IPeng et "al 



( 20061 ) is a little lower than observed for the LWT bursts, 
they only computed a few different cases of hydrogen ig- 
nition. In addition, the recurrence time and burst en- 
ergies do match the expectations for hydrogen ignition 
close to the stability boundary quite well. An interest- 
ing point about this interpretation is that t he sedimen- 



tatio n of heavy elements plays a crucial role (jPeng et al 



Without sedimentation, the CNO abundance re- 
mains close to the solar value at the ignition depth, and 
the energy released by the initial proton captures is small 
giving a very weak flash. With sedimentation, the CNO 
abundance is enhanced by a factor of several, leading to 
more energy release during the initial runaway, giving a 
peak temperature of « 3 x 10 8 K and a peak luminosity 
observable above the accretion luminosity. 

In addition, the interpretation of the LWT bursts as 
hydrogen flashes also provides two explanations for the 
bri ght energeti c burst (£b = 3.6 x 10" 7 erg cm" 2 ) reported 
bv lWolff et all (|2005h from EXO 0748-676 at a low persis- 



tent flux (/ p = 4.8 x 10" 10 erg cm" 2 s" 1 ). First, if the 
hydrogen flash fails to ignite helium, a pure helium layer 
builds up beneath the hydrogen burning shell which will 
eventually ignite by triple alpha reactions. Alternatively, 
since the burst properties indicate that the hydrogen igni- 
tion is occurring close to the transition between unstable 
and stable hydrogen burning, a slight increase in accretion 
rate would lead to stable hydrogen burning and a similar 
accumulation of a pure helium layer. Second, if the ac- 
cretion rate drops, the hydrogen ignition column depth 
becomes large enough t o trigger hel ium ignition resulting 
m a mixed H/He flash. IWolff et alJ (120051) derive a burst 
fhicnce of 3.6 x 10~ 7 erg cm" 2 which translates to a burst 
energy of (1.5-2.5) x 10 39 ergs (depending on whether the 
peak luminosity corresponds to the Eddington limit for 
solar composition or for pure helium; the corresponding 
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distances are 5.9 and 7.7 kpc). Assuming Q nuc = 1.6 MeV, 
appropriate for helium burning to iron group, the implied 
thickness of the fuel layer is « 10 8 g cm~ 2 . A pure he- 
lium layer accumulating beneath the hydrogen shell would 
typically reach much greater thi cknesses before ig niting. 
Therefore, it seems likely that the lWolff et al.l (|2005h burst 
is another example of a hydrogen ignition, but this time 
also igniting the helium. 



6.2. Comparison with theory: SWT bursts 

We have spent some time discussing the origins of the 
LWT bursts in the hope that this might give a clue to 
the origin of the SWT bursts. These have so far evaded 
theoretical explanation, altho ugh possibilitie s have been 
put forward in the literature. iFuiimoto et al.l (|1987l ) sug- 
gested that unburned fuel left over from the preceding 
burst is mixed downwards b y hydrodynamic in s tabili- 
ties driven by rotational shear. IWallace &: Wooslevl (|1984h 
suggested that fresh fuel could be mixed downwards by 
Rayleigh- Taylor instabilities which set in as the layer 
cools. However, neither of these explanations offers a natu- 
ral explanation for the ten minute delay time. This point is 
emphasized dramatically by the discovery of burst triplets 
discussed in this paper. The fact that the same character- 
istic timescale sets the delay between the second and third 
bursts in a triplet as well as the first and second points to 
some physical process with this timescale that can recur 
more than once in succession. 

Mixed H/He bursts at accretion rates close to 0.1 of the 
Eddington rate have been modeled extensively in spherical 
symmetry with large nuclear reaction networks (Woosley 
et al. 2004). However, there has been little work on hy- 
drogen triggered bursts, and so the suggestion that the 
LWT bursts are hydrogen triggered opens up the possi- 
bility that the ten minute phenomenon is connected with 
unstable hydrogen burning. For example, the half-life of 
13 N in the CNO cycle is 9.97 minutes, very close to the ob- 
served timescale. Perhaps as the layer cools following the 
initial flash, an instability driven by proton captures on 
13 C or 14 N occurs once there has been time for enough seed 
nuclei to be produced by the beta-decay of 13 N. However, 
it is not clear how the fact that the timescale between 
SWT bursts changes (see Sect. [5]), possibly with accretion 
rate, would be accommodated in this picture. Further the- 
oretical studies of hydrogen triggered bursts are needed to 
explore the possibility that the SWT bursting has a nu- 
clear physics explanation. 

Even without a good picture of the physics of the SWT 
bursts, we can explore some possibilities by considering 
the idea that the second and third bursts are caused by 
ignition of leftover fuel from the first burst. One way that 
this might happen is that the ignition conditions for the 
first burst are always the same, but that sometimes incom- 
plete burning occurs, leaving behind unburned fuel that 
later reignites leading to double or triple events. However, 
this is inconsistent with the observation that, at a given 



wait time, the fluence is the same for a Dl, a Tl or a SI 
burst (the LWT bursts have all the same alpha values). In 
addition, the fact that the wait time to the next event is 
longer after a double or triple (regardless of the type of the 
next event) also suggests a different picture. Imagine that 
after the first burst (either SI, Dl, or Tl), a fraction / of 
the fuel layer is left unburned. Then, in double or triple 
events, this leftover fuel layer is burned after a delay time 
of « 10 minutes, whereas in single events, the unburned 
fuel survives, and leads to early ignition of the next event. 
The difference between singles and doubles/triples here 
is not whether there is any unburned fuel, but instead 
whether the unburned fuel is able to ignite on the ten 
minute timescale. 

One way to test this idea is to look for consistency be- 
tween the wait times and the burst fluences. The wait time 
after a double (triple) is « 37% (55%) longer than the wait 
time after a single, suggesting that a fraction f a » 37% 
(55%) of the available accreted fuel is not burned dur- 
ing the first burst. If this unburned fuel was of the same 
composition as the rest of the layer and burned subse- 
quently in the second/third bursts, we would expect the 
double/triple events to have a larger fluence by the same 
factor f a . Now, the double (triple) events have larger flu- 
ences than the singles, but only by a factor /(, of w 17% 
(25%). This indicates that the layer that burns in the sec- 
ond/third bursts has a lower energy per gram, for exam- 
ple a smaller hydrogen content than the accreted fuel that 
burns in the first burst. This suggests that in fact the "un- 
burned" fuel layer undergoes some hydrogen burning into 
helium during the first burst. Interestingly, the ratio f a /fb 
is « 2.2 for both the doubles and the triples and is < 2.5 
which is the value of the ratio between the energy released 
by an hydrogen nucleon burning to the iron group (4 MeV) 
and the energy released by an helium nucleon burning to 
the iron group (1.6 MeV). The burst lightcurves also sup- 
port the idea that the second/third bursts have a smaller 
hydrogen fraction since they have a much shorter tail than 
the first bursts. 

These properties would be naturally explained if the 
first burst in the sequence involved unstable hydrogen 
burning into helium. If hydrogen ignition is occurring near 
the transition from unstable to stable hydrogen burning, 
then it is not clear whether helium is burning at all in 
the first flash. Certainly, it is not able to ignite unsta- 
bly at low column depths but maybe it burns slightly at 
the peak temperature. The first flash would thus leave a 
large part (if not all) of the available helium unburned and 
further produce some. A subsequent second/third burst 
would take place in an helium-rich environment. The fact 
that little material can be accreted within ten minutes 
suggests that the SWT bursts are also hydrogen triggered 
since lower column depths are required than for helium ig- 
nition. More theoretical work on hydrogen ignition bursts 
is needed to test these ideas. 

One important question concerns the occurrence rate 
of quadruple events. We did not observe a qu adruple 
burst from EXO 0748-676. iGallowav et~ai1 (|2006l ) report 
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one quadruple event from 4U 1636-536 in the RXTE burst 
catalog. The occurrence rate of quadruple versus triple 
events (as well as triple versus double events) could pro- 
vide an important constraint on theoretical models. 

Conclusion 

This is the first time that unambiguously and repeat- 
edly triple bursts are detected in an accreting neutron 
star, EXO 0748-676. The source was already known to fre- 
quently exhibit double bursts, where two bursts are sep- 
arated by only 12 minutes. Now the same time scale is 
seen in triple bursts. This recurrence time is too short to 
accrete a sufficient amount of fuel, which indicates that 
there must have been fuel left unburned in the previous 
burst. We suggest that after each burst there is an amount 
of fuel left unburned, which may or may not reignite after 
12 minutes, in the former case producing a double or triple 
burst. We find evidence that this unburned fuel yields less 
energy per gram, pointing to a lower hydrogen content 
than for the fuel of the primary bursts. 

The recurrence times of 3 hours and the low alpha val- 
ues are consistent with mixed hydrogen/helium burning 
at 0.1 Eddington, as seen for example in GS 1826-24. The 
persistent luminosity of EXO 0748-676 is consistent with 
this accretion rate if the X-ray emission from the system 
is anisotropic, as predicted for large inclination angles. 
However, the expected energies and peak luminosities are 
somewhat larger than observed in EXO 0748-676. In ad- 
dition, this burning regime has been well-studied theoret- 
ically and does not predict the ten minute bursting phe- 
nomenon. New physics, such as mixing of fuel to deeper 
layers, is required to explain the multiple events within 
the mixed hydrogen/helium burning regime. 

Hydrogen ignition at a lower accretion rate may pro- 
vide a more natural explanation for the burst behavior 
of EXO 0748-676. Hydrogen ignition bursts are not well- 
studied, but the properties of the bursts from EXO 0748- 
676 ag ree quite well with the recent work by IPeng et al. 
( 20061 ). In particular, the low peak luminosity and burst 
energies of the LWT bursts, and the low persistent lu- 
minosity are consistent with burning in this regime. This 
opens up the exciting possibility that the ten minute re- 
currence time bursts are in fact a natural consequence of 
hydrogen ignition bursts. 
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Furthermore, since an eclipse lasts ~500 s, which is less 
than the typical separation of the bursts in a multiple 
burst event, it is not possible either to have missed a com- 
plete doublet or triplet. From the inspection of Fig. IA.1I 
we can therefore conclude that it is unlikely that any SI, 
Dl, Tl, D2, T2 or T3 burst has been missed because of 
an eclipse. 

However, bursts occurring close to (less than 25 min) 
the beginning or the end of an observation, or close to a 
short instrumental gap such as occasionnally present in 
the data, cannot have their type determined certainly be- 
cause another burst of the same event could have occurred 
during the gap and been missed. In such cases, the differ- 
ent possible types are indicated in Table [A~T1 starting with 
the one attributed to the burst in this analysis. 



Appendix A: Light curves & bursts lists 

Fig.EZDshows the EPIC PN light curves of EXO 0748-676 
with a binning of 60 s during XMM-Newton revolutions 
692 to 719. Panels a and b show the 5-10 keV and the 
0.3-5 keV energy band, respectively. Panel c shows the 
hardness ratio. Times are not barycentre-corrected. A cy- 
cle number tracing the eclipse timing is indicated on the 
top axis. This was determined from a reference time of 
13.79 h on 2003 September 19 (or XMM-Newton time of 
1.8036643e + 08 s) estimated as the mid-time of the first 
eclipse observed during revolution 692, a nd using an or- 
bital period of 3.824 hr (|Wolff et al.ll2002h . Integer values 
of the cycle number correspond to phase 0, or estimated 
eclipse mid-times. 

The peak time, wait time and type of all the bursts are 
listed in Table [A~T1 Since the wait time of the first burst 
of each observation in the series cannot be determined, 
a lower limit is given as the separation between the burst 
peak time and the start time of the observation. An upper 
limit is given as the separation between the burst peak 
time and peak time of the last burst observed during the 
previous observation of the series. 

We note that because of scattering, residual emission 
is present during eclipses at a level of ~4% of the source's 
persistent flux. This explains why an X-ra y burst was once 
detec ted during an eclipse with EXO SAT ( Gottwald et alj 
19861 ). This implies that any SI, Dl or Tl (ie. strong) 



burst occurring during an eclipse of the XMM-Newton 
observation would have been detected. Since we do not 
detect any burst during eclipses, it appears unlikely that 
we have missed an SI, Dl or Tl burst because of eclipses. 
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Fig. A.l. EXO 0748-676 as observed by XMM-Newton EPIC PN during revolutions 692 (top) and 693 (bottom), a) 
5-10 keV light curve. The bursts are numbered, b) 0.3-5 keV light curve, c) Color (counts in the 5-10 keV band 
divided by counts in the 0.3-5 keV band) as a function of time. The binning time is 60 s in each panel. The cycle 
number is indicated on the top axis, with the estimated mid-eclipse time of the first eclipse observed during revolution 
692 taken as cycle number 1. 




Fig. A.l. Continued for XMM-Newton revolutions 694 (top) and 695 (bottom). The dashed line in panel c of the 
bottom plot indicates the time of a burst (numbered 39) detected by RGS just before the start of the EPIC PN 
observation. 
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Table A.l. Bursts timing. The columns indicate the burst number, the XMM-Newton revolution number at which 
the burst occurred, the burst peak time determined from the EPIC PN 5-10 keV 1 s resolution light curve (except for 
burst 39 detected only by RGS), the wait time and the burst type. S means single burst. Dl and D2 mean first and 
second burst in a doublet. Tl, T2 and T3 mean first, second and third burst in a triplet, respectively. When the burst 
type cannot be determined certainly, we indicate the different possible types, starting with the one attributed to the 
burst in this analysis. 



Burst Rev. Peak time (UT) t wait (h) Type 

1 692 Sep 19 15:13:58 1.62 < SI 

2 ' 19:52:12 4.64 SI 

3 23:04:31 3.21 SI 

4 20 01:36:43 2.54 Dl (T2) 

5 01:46:28 0.162 D2 (T3) 

6 05:21:42 3.59 SI (D2) 

7 07:58:35 2.61 Dl 

8 08:10:17 0.195 D2 

9 11:31:58 3.36 SI 

10 13:57:49 2.43 SI (Dl, Tl) 

11 693 21 13:53:56 0.259 < < 23.9 SI (D2, T3) 

12 16:08:03 2.24 SI 

13 18:32:55 2.41 Tl 

14 18:43:22 0.174 T2 

15 18:55:05 0.195 T3 

16 23:15:33 4.34 SI (D2) 

17 22 02:09:52 2.91 SI 

18 04:26:17 2.27 Dl 

19 04:35:11 0.148 D2 

20 07:10:07 2.58 SI (D2) 

21 09:23:54 2.23 SI 

22 11:37:24 2.23 Tl 

23 11:51:07 0.229 T2 

24 12:04:08 0.217 T3 

25 694 23 12:32:29 1.83 < < 24.5 SI 

26 15:05:56 2.56 SI 

27 17:25:09 2.32 SI (D2) 

28 20:00:08 2.58 SI 

29 22:34:36 2.57 Dl 

30 22:45:51 0.187 D2 

31 24 02:10:45 3.41 Dl 

32 02:26:15 0.259 D2 

33 06:03:06 3.61 SI 

34 08:50:34 2.79 Tl 

35 09:05:55 0.256 T2 

36 09:16:31 0.177 T3 

37 13:36:43 4.34 SI 

38 16:37:55 3.02 SI (Dl, Tl) 
a 39 695 25 17:28:20 0.036 < < 24.8 Dl (T2) 

40 17:41:31 0.220 D2 (T3) 

41 21:34:54 3.89 SI 

42 26 00:28:00 2.89 Dl 

43 00:46:20 0.306 D2 

44 04:28:35 3.70 Dl 

45 04:44:35 0.267 D2 

46 08:51:35 4.12 SI 

47 11:41:43 2.84 SI 

48 708 Oct 21 11:51:35 1.82 < < 600 SI (Dl) 

49 14:12:10 2.34 SI (Dl) 

50 16:30:24 2.30 Tl 

51 16:39:49 0.157 T2 

52 16:51:05 0.188 T3 

53 19:59:49 3.15 Dl 

54 20:18:52 0.317 D2 

55 22 00:28:46 4.17 SI 

56 710 25 20:58:42 1.65 < < 92.5 SI 

57 23:48:23 2.83 Dl 

58 26 00:01:37 0.221 D2 

59 03:11:08 3.16 SI (Dl) 

60 05:38:40 2.46 SI (D2) 

61 07:41:16 2.04 Dl 

62 07:53:28 0.203 D2 

63 10:58:02 3.08 Dl 

64 11:13:05 0.251 D2 

65 719 Nov 12 08:35:44 0.191 < < 405 SI (D2, T3) 

66 11:42:47 3.12 SI 

67 15:23:00 3.67 SI 

68 18:31:01 3.13 SI 

69 21:17:46 2.78 Dl 

70 21:26:09 0.140 D2 

71 13 01:04:58 3.65 SI 

72 03:39:57 2.58 Tl 

73 03:51:47 0.197 T2 

74 04:03:36 0.197 T3 

75 07:27:43 3.40 Dl 

76 07:39:37 0.198 D2 



a detected by RGS only, before EPIC cameras were turned on. 



